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Abstract 

In this work, we derive coarse-grained finite-temperature theory for a Bose condensate in a one- 
dimensional optical lattice, in addition to a confining harmonic trap potential. To construct the theory for the 
condensate and noncondensate in a periodic lattice potential, the difficulty arises due to the rapid variation 
in the position by a lattice potential, compared to the length scale of the harmonic potential. In order to 
overcome this difficulty, we need some coarse-graining procedure for the lattice potential. We start from a 
two-particle irreducible (2PI) effective action on the Schwinger-Keldysh closed-time contour path. In prin- 
ciple, this action involves all information of equilibrium and non-equilibrium properties of the condensate 
and noncondensate atoms. By assuming the ansatz for the variational function, i.e., the condensate order 
parameter in an effective action, we derive a coarse-grained effective action, which describes the dynam- 
ics on the length scale much longer than a lattice constant. Using the variational principle, coarse-grained 
equations of motion for the condensate variables are obtained. These equations include a dissipative term 
due to collisions between condensate and noncondensate atoms, as well as noncondensate mean-field. As 
a result of a coarse-graining procedure, the effects of a lattice potential are incorporated into equations of 
motion for the condensate by an effective mass, a renormalized coupling constant, and an umklapp scatter- 
ing process. To illustrate the usefulness of our formalism, we discuss a Landau instability of the condensate 
in optical lattices by using the coarse-grained generalized Gross-Pitaevskii hydrodynamics. We found that 
the collisional damping rate due to collisions between the condensate and noncondensate atoms changes 
sign when the condensate velocity exceeds a renormalized sound velocity, leading to a Landau instability 
consistent with the Landau criterion. Our results in this work give an insight into the microscopic origin of 
the Landau instability. 



1 



I. INTRODUCTION 



Recent extensive researches on ultracold atomic gases in optical lattices by experimental and 
theoretical approaches have revealed the nontrivial nature of many-body quantum systems, [1,2] 
such as the superfluid-Mott insulator transition, [3, 4] stability of superfluidity, [5-9] and the 
Josephson effect. [10] Because of the ease of fine-tuning of experimental parameters, optical 
lattices are also used as testing grounds for many-body theory, [11] including non-equilibrium 
phenomena, [12-17] which are usually very difficult to analyze in traditional solid state systems. 

Among a number of startling behaviors of the Bose condensate, superfluidity is one of the most 
fascinating phenomena. A long time ago, Landau showed that the superfluid state is stable as long 
as a velocity of a superfluid is smaller than a critical velocity, above which elementary excitations 
are spontaneously produced, making the superfluid state unstable. [18, 19] His argument relied 
only on the energy and momentum conservation and the Galilei transformation. When a conden- 
sate is set in an optical lattice potential, the breakdown of superfluidity becomes more complicated 
due to the competition between the inter-atomic interaction and the periodic lattice potential, as 
observed experimentally. [5, 7, 9, 20] Theoretically, two different types of instability have been 
discussed within the Gross-Pitaevskii (GP) equation with a periodic lattice potential. [21-28] One 
is the Landau (or energetic) instability, which occurs when the excitation energy becomes neg- 
ative. It is this instability that Landau originally argued. The other is the dynamical instability, 
which occurs when the excitation energy possesses the imaginary part. In this type of instability, 
the effect of the lattice potential is to couple an (unphysical) antiphonon to a phonon by the first 
order Bragg scattering, leading to the dynamical instability. [21, 22, 26] It is important to note that 
the dynamical instability exists at zero temperature, and thus can be understood within the usual 
zero-temperature GP theory. In contrast, however, the Landau instability has been found to occur 
at finite temperatures. [9] As compared with the dynamical instability, the essential role of the 
lattice potential in the Landau instability is to pin the incoherent thermally excited noncondensate, 
while the condensate can coherently tunnel through the lattice potential. [20] Thus the thermally 
excited atoms trapped by the optical lattice play a role of obstacles to the condensate, giving rise 
to dissipative effects. For the above reason, one cannot study microscopic mechanisms of the Lan- 
dau instability by using the zero-temperature GP equation. One should include the effect of the 
thermal cloud into the theory. 

In order to discuss the Landau instability in the periodic lattice, one needs a finite-temperature 
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microscopic theory for the Bose-condensed gas including effects of the lattice potential. As far as 
very low temperature regime is concerned, the GP equation has succeeded in describing a trapped 
Bose-Einstein condensate. [29-31] At finite temperatures, however, the presence of the noncon- 
densate in addition to the condensate makes the GP description insufficient. In order to incorporate 
the dynamics of the noncondensate, a number of papers derived generalized GP equations at finite 
temperatures, which includes effects of noncondensates by mean-field and collisional exchange 
between condensate and noncondensate atoms, and a quantum kinetic equation for the noncon- 
densate. [32-36] A resulting two-component many-body system leads to non-trivial phenomena 
such as nucleation and evaporation of condensates, [34-38] and damping of collective modes, [39- 
42] which do not appear in the GP theory. It is thus natural to anticipate that incorporation of the 
periodic lattice potential into the coupled many-body system of the condensate and noncondensate 
atoms will lead to interesting new physics. The main purpose of the present paper is to show one 
of the crucial effects due to thermally excited noncondensate atoms in optical lattices by focusing 
on the breakdown of superfluidity. 

In the present paper, for investigating non-equilibrium dynamics of such a two-component 
many-body system in the periodic optical lattice potential, we construct a finite-temperature the- 
ory for the condensate in the one-dimensional optical lattice, in addition to a confining three- 
dimensinal trap potential. For this purpose, we start from the 2PI effective action [43] with the 
Schwinger-Keldysh closed-time path formalism. [44-46] There are two advantages to use a func- 
tional integral formalism. The one is that one can derive an action in the non-equilibrium quantum 
field theory by controllable approximations. [13, 14, 47, 48] The other is that one can introduce an 
ansatz as the variational function, i.e., the order parameter of the condensate in the effective ac- 
tion and perform coarse-grained approximation in a well-defined manner since the action involves 
integrations over position. Making coarse-grained procedure, one can include the effects of the 
lattice potential effectively. By means of the coarse-grained effective action, we obtain coarse- 
grained equations of motion for the condensate variables at finite temperatures including the effect 
of the optical lattice. As an application of our finite-temperature coarse-grained theory, we inves- 
tigate the microscopic mechanism of the breakdown of the superfluidity having experiments such 
as reported in Ref. 9 in mind. 

Recent papers 49-5 1 reported the earlier attempts to study the breakdown of superfluidity in an 
optical lattice, focusing on microscopic mechanisms. When the condensate has finite velocity, the 
damping rate of the collective oscillation, which stems from the collisional or Landau damping 
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processes, can change its sign at a critical velocity. This means that the inverse process of damp- 
ing occurs at the critical velocity, resulting in spontaneous production of elementary excitations. 
Moreover, the increase of excitations in time as a result of inverse damping process suggests the 
breakdown of the stable superfluid state. In Refs. 49-51, starting from the one-dimentional Bose- 
Habbard model, the authors derived equations of motion for the condensate order parameter at 
finite temperatures, which include the effect of the noncondensate atoms through mean-field inter- 
actions as well as collisions. From these equations of motion, the authors calculated damping rate 
of collective modes (phonons) in the collisionless regime of interest and showed that the damping 
rate changes its sign at a critical velocity. This instability is shown to coincide with the Landau 
instability. 

In contrast to the previous works, [49-51] in the present paper we use a coarse-grained the- 
ory. It will be shown that the coarse-grained formalism developed here describes the breakdown 
mechanism in a more transparent way. From the condition of the negative damping rate, we auto- 
matically obtain both the negative excitation energy and the Landau criterion, which is modified 
by the lattice potential. 

The present paper is organized as follows. In Sec.II, we briefly review the non-equilibrium 
quantum field theory that consists of the 2PI effective action with the Sch winger- Keldysh closed- 
time formalism. In Sec. Ill, for illustration of the approximation, we derive a generalized GP 
equation for the condensate, which can be written as the hydrodynamic equations in terms of the 
local condensate density and the superfluid velocity. In Sec. IV, we derive the coarse-grained ef- 
fective action for the Bose gas in optical lattices in terms of coarse-grained macroscopic variables. 
Using the variational principle, we obtain coarse-grained equations of motion of the condensate 
variables suitable for describing the condensate at finite temperatures in an optical lattice. In 
Sec. V, in order to show the usefulness of our theory, we apply coarse-grained generalized GP 
hydrodynamic equations to discuss the breakdown of the superfluidity and give the microscopic 
origin of the Landau instability. 

H. 2PI EFFECTIVE ACTION AND NON-EQUILIBRIUM QUANTUM FIELD THEORY 

We consider a Bose-condensed gas trapped in a one-dimensional optical lattice in addition to 
the harmonic potential: 

Kxt(r) = \4o(r) + KptW, (1) 
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KptO) = sE R cos 2 , (3) 

where u x , to y , lo z are frequencies of the harmonic trap potential, s is the dimensionless parameter 
describing strength of the lattice potential, E R = h 2 /2md 2 is the recoil energy, and d is the lattice 
constant. A Bose gas in the external potential represented by Eq. (3) is described by the following 
Lagrangian density 

-Kxt(r)|^(r,t)r-||^(r,t)| 4 , (4) 

where ij) is the scalar field describing Bose atoms. We have assumed that the interaction between 
atoms is represented as a short-range peudopotential with the coupling constant g, which is related 
to the s-wave scattering length of atoms through g = 4irh 2 a/m. 



A. Generating Functional in the Non-equilibrium Quantum Field Theory 

An efficient way to treat non-equilibrium dynamics [47] is using the two-particle irreducible 
(2PI) effective action [43] with the Sch winger- Keldysh closed-time-path formalism. [44-46] This 
formalism provides a powerful starting point for systematic approximations in the non-equilibrium 
quantum field theory. [13, 47, 48] For the ultracold atomic Bose gases, Rey et al. first applied 
this method to study various approximations and to develop the quantum kinetic theory for the 
condensate in optical lattices. [13, 14] From the 2PI effective action derived below, one can obtain 
a generalized GP equation for the condensate order parameter including effects of noncondensate 
atoms as a dissipative term and the noncondensate mean-field, and some kind of a kinetic equation 
for noncondensate atoms on an equal footing. In this section, we briefly review this technique (see 
Ref. 47 for more details). We use units where % = 1 in this section. 

In the quantum field theory, any correlation functions, which contain all information about a 
non-equilibrium many-body system, are obtained from a generating functional. As long as the 
initial density matrix is approximated by the Gaussian form, the generating functional in the non- 
equilibrium field theory can be written as a functional integral [47, 48] 



Z[J,K] = Jv^exp i(s[4>] + J^il> + -^hKil> 



(5) 
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FIG. 1: The Schwinger-Keldysh contour path 



where we have introduced matrix notation for the scalar field ip and the source field J 

i>(r,t) = ^(r,t),r(r,t)Y = (Mr,t),Mr,t)Y, 

J(r,t) ee (J(r,t),J*(r ! t))*EE(J 1 (r,t),J 2 (r,t)) t , (6) 

and suppressed the integration over space and time variables. A 2 x 2 matrix nonlocal source field 
K is added to the action in order to obtain equations of motion for two-point correlation functions 
(Green's functions). The classical action is defined by the Lagrangian density (4) 

S[ip] = J dr J dtC(r,t), (7) 

where the subscript C of the time-integration means that the integral is performed on the 
Schwinger-Keldysh contour path, which extends from the initial time t to the finite time t > t , 
and back from t to to (Fig.l). From the generating functional (5), one can define a more useful 
generating functional for the connected Green's function 

W[J,K} = -i\nZ[J,K}. (8) 

The derivative of the generating functional W with respect to the source field J gives the conden- 
sate order parameter 

™ = *<(M) « = 1,2). (9) 
dJi{r,t) 

This order parameter $(r, t) is the average of the original complex field -0 taken by the action (7) 



$(r,t) ee ($(r,t),$*(r,t))* 

EE (^(r, t),$ 2 (r,t))' 

ee / Dip tp exp 

ee (VM)>. (10) 



On the other hand, the derivative of W [J, K] with respect to the nonlocal source field K gives 

5W[J,K] 1 



6K(r,t; r',t') 



[H(r,t;r',t')+iG(r,t-r',t')} 



(ID 



Here, we have defined the condensate matrix Green's function, H(r, t; r', t'), and the nonconden- 
sate connected matrix Green's function, G(r, t; r', t')\ 



H(r,t-r',t' 



$(r,t)$*(r',t') $(r,t)$(r',t') 
$*(r,t)$*(r',t') $*(r,i)$(r',t') 



(12) 



iG(r,t; r',t" 



(13) 



In the noncondensate Green's function, giving by (13), we introduced the noncondensate field: 



= (^(r,t)-$(r,t),^(r,t)-$*(r,t))' 



(14) 



We note that averages over the fields are automatically time ordered on the closed-time path in the 
functional integral. 

We now define the effective action, which is the generating functional for the two-particle 
irreducible vertex functions, through the Legendre transform 



r[*, G] = W[J, K] - J f $ - -&K3> - -Tr [GK] . 
Following Ref. 43, one obtains the expression for the effective action r[<fr, G] as 

r[*,G] = + ^Tr In G 1 + ^Tr [D^G] +T 2 [^,G} + Const., 



where D 1 is the classical inverse propagator matrix defined by 

5 2 S[^} 



(15) 



(16) 



(17) 



with the 2x2 matrix elements 

£>ii 1 (r,i;r' J 



5 2 S[<$>] 



5$*(r,t)5$(r',t') 
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5(r - r')6(t - t') 



D^(v,t;r',t') 



D^(r,t; r',t>) 



D^(r,t;r',t') 



-2^|$(r)| 2 (5(r-r'), 

5 2 S[<f>] 
5$(r,t)5$(r',t') 

-2g [$(r)$(r')] 2 <5(r-r'), 

5 2 S[<f>] 
6$*(r,t)5$*(r',t') 
-2g [$*(r)$*(r / )] 2 5(r-r / ) 

5 2 S[$] 
5$(r,t)5$*(r',t') 
. d 



(18) 



(19) 



(20) 



<5(r - r')5(t - t') 



-2#|$(r)| 2 5(r-r'). 
In Eqs. (18)— (21), Hq(t) is defined by a one-body part 

ft (r) = --^V 2 + Kx t (r). 
2m 



(21) 



(22) 



The trace and logarithm in Eq. (16) is defined by the functional integral. T 2 [&, G] in Eq. (16) 
consists of two-particle irreducible vacuum diagrams (the diagrams that cannot be disconnected 
by cutting two propagator lines) with full propagators G. The vertices are determined by the 
interaction term Si n t[<fr, i/>], which is the part higher than second order in ip of the expansion for 
the action S[$> + ■»/>]. 

In this paper, we consider a relatively high-temperature regime by treating noncondensate 
atoms within the Hartree-Fock approximation, neglecting off-diagonal components of the Green's 
functions. The resulting effective action is given by 



r[*,G] = r*[*,G] + r G [*,G], 



(23) 



where 



S GGP [<S>,G] 

J dr jf dt$*(r,f){i^ - H (r) - | |$(r,t)| : 
G n (r,t;r,t) + G 22 (r,t; r,t) 



r„r*,Gi 



- w 
i 



(24) 







j dr J dt 


4-«o(r) 





iG u (r,t;r,t) 
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+- / 



4 -«•(') 



iG 22 (r,t; r,t) 



+^TrlnG~ 1 . (25) 

In Eq. (23), the matrix Green's function contains only diagonal components because we will use 
the Hartree-Fock approximation for the noncondensate. The subscript GGP of Sggp m Eq. (24) 
indicates that Sggp will be shown to lead to the generalized GP equation. 

Taking the derivative of the effective action (23) with respect to $ and G, one obtains ST / 5<& = 
—J — K<& and 5T/5G = —iK/2, respectively. In a real physical state, the artificial external 
fields J and K should vanish. This requirement yields equations of motion 5T/5&*(r, t) = 
5S GGP /5$*(r,t) = 0: 



zJU(r,t) =W„(r)$(r,t)+0 



®(r,t)\ 2 + iGu(r,t;r,t) 



+ iG 22 (r,t; r,t) 

and 5T/5G u (r,t; r',t') = 0: 

G n X (r , *; r', = D n l (r , t; r', t') - E n (r, i; r', if) . (27) 

Eq. (26) is the equation of motion for the condensate order parameter, while Eq. (27) is the non- 
equilibrium Dyson equation for the noncondensate atoms with the proper self-energy defined by 



S n (r,t;^t0^2z (28) 



Note that since the time integration in Eq. (27) is defined on the Schwinger-Keldysh contour path, 
this equation is generalization of the usual Dyson equation to the non-equilibrium systems. After 
projecting the time-integration on the Schwinger-Keldysh contour path onto the real-time contour, 
Eq. (27) leads to a kinetic equation for the noncondensate distribution function. [52] 

In order to obtain the 2PI part T 2 , one has to perform approximations suitable for the physical 
problem under consideration by truncating diagrammatic expansions. We show diagrams for T 2 
used in the present analysis in Fig. 2, where two- and three-loop vacuum diagrams are illustrated. 
The vertices are specified by the interaction Si n t[3?, Analytical expressions corresponding to 
Fig. 2 are 

r 2 [<&, G] = if >[*, G] + T 2 2:c) [$, G] + if nc) [$, G], (29) 
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where 



Gu(r,t; r,t)Gn(r,t; r,t) 
+ 6Gn(r,t; r,t)G 22 (r,t;r,t) 
+ G 22 (r,t;r,t)G 22 (r,t;r,t) , 



drdr' / dtdt' 



4 

xG 22 (r,t; ^,0^22^, t; r',t' 
+ # 22 (r,t; r^G^r,*; r',t' 



iZnfotjr^Gnfot; r',t' 



,(2:nc) 



drdr' / dtdt' 



G 11 (r,t; r',t')G u (r,t ]r ',t' 



x G 22 (r,t; r',t')G 22 {r,t;r',t' 



From these equations, self-energies are derived by using the relation (28): 



Sff(r,t;r',t') 



ig 

2 



Gii(r, t:r', t') + G 22 (r,t: r't 



5(t - t')6(r - r') 



+ igG 22 (r,t ]r ',t')5(t-t')6(r- r '), 



ESf)(r,t;r',t') = -|- 



#u (r', f ; r, t)G n (r, t; r', f)Gn (r, t; r', f) 
+ 2ff u (r, *; r\ f')Gn(r, *; r', 0^11^, f; r, t) , 
Eg :nc) (r, *; r', *') = -2ig 2 G^T, *; r', t')G u (r, *; r', t^G^r' , f; r, t). 



(30) 



(3D 



(32) 



(33) 

(34) 
(35) 



Before closing this section, we mention some approximations for the effective action. [13, 14, 
47, 48] The GP equation is obtained by retaining only the classical action S[<&] in Eq. (16). This 
corresponds to the mean-field approximation appropriate at zero temperature. When one retains 
all terms except T 2 , this yields the Bogoliubov or one-loop approximation. Including T 2 up to first 
order in the coupling constant g, one obtains the time-dependent Hartree-Fock-Bogoliubov equa- 
tions. [32] One needs to consider higher-order approximation to include the multiple scattering 
effect. [13, 47, 48] It will be shown that the multiple scattering gives rise to the dissipative term in 
the equation of motion for the condensate and the colllision integrals in the kinetic equation. 
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FIG. 2: 2PI diagrams up to second order in the coupling constant; the double-bubble (a), the setting sun (b), 
and the basketball (c). 

III. MICROSCOPIC FINITE-TEMPERATURE THEORY FOR THE BOSE GASES: A REVIEW 

It is instructive to derive the generalized GP equation and hydrodynamic equations for the con- 
densate at finite temperatures, which have been derived in Ref. 35, within the 2PI formalism. In 
this section, we give a derivation of a generalized GP equation by specifying diagrams to approx- 
imate T 2 . In Sec. IV, we will use the technique discussed in this section to derive generalized 
hydrodynamic equations including the effect of optical lattices. 

A. Generalized Gross-Pitaevskii Equation 



Using Eqs. (30), (31), and (32), one obtains the analytical expression of the last term of Eq. (26) 



as 



dr' / dt 



ST 2 [<f>,G] = _g^_ 
6$*{r,t) 2 

x Gn(r, *; r', t')G n (r : t; r', t')G n (r', t'; r, t)$(r', t') 



dr' J dt'F(r,t; r', t')<5>(r' ,t') 



(36) 



where we have used the relation G 2 2(r, t; r', t') = Gu(r', t'; r, t). For later convenience, we have 
introduced the following quantity 



F(r, t- r', t') = 9 -G{v, t- r', t')G(r, t; r', t')G{r', t'; r, t). 



(37) 



Here and hereafter, we use the notation G(r, t; r', t') instead of Gii(r, t; r', t' ). In Eq. (36), it 
should be noted that the time integration is defined on the Sch winger- Keldysh contour path. In 
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order to perform the integration in Eq. (36) explicitly, one has to project the time integration 
on the Schwinger-Keldysh contour path onto the real time axis. By virtue of the principle of 
causality, the integrand in Eq. (36) is replaced with the retarded counterpart. After using the 
Langreth theorem, [52] one obtains 



- / dt'F(r,t;r',t')<l>(r',t f ) 

J c 



9 2 



- f°° rft , b (+) (r,t;r / ,t , )G' < (r,t;r , ,t')G' < (r / ,t / ;r,t) 

J — oo 
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+ G< (r, t- r', t')G {+) (r, t; r', f)G< (r', t>; r, t) 
+ G {+) (r, t; r', t')G M (r, t; r', t')G< (r', 1'; r, t) 
+ G<(r, t; r', t)G K {r, t; r', t')G^\r', t'; r, t) 
x$(r',t'), (38) 
where we have introduced the lesser, greater, retarded, and advanced Green's functions as 

G < (r,t;r',t') = -i(r(r>,t>)4>(r,t)), (39) 

G>(r,t;r',t') = -i$(r, ^*(r', 0), (40) 

G (+) (r, t; r, = - i')([^(r, t), ^*(r\ *')]), (41) 

(?H ( r , t ; r\ t') = i0(f - t) ($(r, t), ^(r', *')]>• (42) 

Eq. (38) involves terms that are nonlocal in space and time, which make it difficult to solve 
the equation. As in Refs. 35, 36, we assume that the macroscopic variables vary slowly in space 
and time. We thus approximate the condensate order parameter near the specific position and time 

(r, t) as 



$(r',t') = ^n c {j\t')e ie{r '^ 

~ y / n c (r,t)e i[9(r ' i)+ate(r ' <)(i '~' )+w(r ' i) ' {r ' _r)1 

= $(r, t)e~^ a; ^ r ' t - )( '*'~ < ' ) ~ kc< ' r ' < ' ) ^ r '~ r ' ) ', (43) 

where n c (r, t) and 9(r,t) are the condensate density and the phase of the order parameter, re- 
spectively. The condensate frequency and wavevector are defined by u c (r,t) = —d t 9(r,t) and 
k c (r, t) = V0(r, t), respectively. Next, we perform the gradient expansion for the noncondensate 
Green's functions in Eq.(38) in order to separate the scale. For this purpose, we introduce the 
relative coordinate and time and the center of mass coordinate and time 

/ t-. r — r ' 
r = r - r , R = , 
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t = t-t, T ^Af- (44) 

Here the relative coordinates, (r, t), describe the microscopic "fast" dynamics and are treated 
exactly, while the center-of-mass coordinates, (R, T), describe macroscopic "slow" dynamics and 
are treated semiclassically. In order to separate out into the variables describing "slow" and "fast" 
processes, we introduce the Wigner representation, which is defined by the Fourier transforms of 
the relative coordinates 

G(v, i; R, T) = / ^ / ^ e^~^G(k, u>; R, T). (45) 

After performing the approximation fro the condensate (43) and the gradient expansion for the 

noncondensate Green's functions through the Wigner transformation (45), Eq (36) becomes 

ST 2 [^,G] 
<J$*(r,t) 

dt'F(r,t: r'A'Wr'A') 



.g 2 f . f dki dk.2 dk 3 r dui\ du)2 cL> 3 
i — / dr 



(2tt) 3 (2vr) 3 (2tt) 3 J 2tt 2tt 2tt 

x e i(ki+k2-k3-k c )-r 

uj c — uji — uj 2 + cl> 3 + ir] 
x I [G > (k 1 , Wl ; R, T) - G < (k 1 , wi; R, T)] G<(k 2 , u 2 ; R, T) 
xG < (k 3 , W3 ;R,T) 
+ G<(k!, Wl ; R, T) [G>(k 2 , ^ 2 ; R, T) - G<(k 2 , o; 2 ; R, T)] 

xG < (k 3 ,^ 3 ;R,r) 
+ [G > (k 1 , u x] R, T) - G < (k 1 ,u l] R, T)] 

x [G>(k 2 , co 2 ; R, T) - G<(k 2 , o; 2 ; R, T)] G<(k 3 , o; 3 ; R, T) 
-G < (k 1 ,^ 1 ;R,T)G < (k 2 ,^ 2 ;R,T) 

x [G > (k 3 ,^ 3 ;R,T)-G < (k 3 ,^ 3 ;R,T)]|$(R,T). (46) 

We now define the non-equilibrium spectral function: 

A(k, uj; R, T) = % [G> (k, uj; R, T) - G< (k, uj; R, T)] 

= -2lmG R (k,uj;R,T). (47) 

When we introduce a new unknown function f(k,u; R, T) by [53] 

zG<(k, cj; R, T) = A(k, uj; R, T)/(k, uj; R, T), (48) 
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iG> (k, lu; R, T) = A(k, w; R, T) [1 + /(k, cj; R, T)} , (49) 

the relation Eq. (47) is satisfied. The non-equilibrium spectral function can be obtained from the 
equation of motion for the retarded Green's function, G R . In the quasi-particle approximation, the 
non-equilibrium spectral function is given by [54] 

A(k, u- R, T) ~ 2tt5(cj - e(k; R, T)/%), (50) 

where the Hartree-Fock spectrum for the noncondensate is defined by 

e(k; R, T) = — + 2^[n c (R, T) + n(R, T)] + K xt (R), (51) 

with the noncondensate density being defined by n(r, i) = (i/>*(r, t)^(r, £)). In the quasi-particle 
approximation (50), the unknown function / is found to be equivalent to the Wigner distribution 
function, which is defined by 



f w (k,R,T) = J ^-iG<(k,u;R,T), (52) 



where 



iG < (k,w;R J T) 

= |rf te -^|dre 4k - r (^(R+^T+0^(R-^T-0). (53) 

The Wigner function is a quantum counterpart of the classical phase-space distribution function. 
Using Eqs. (48)~(51) in Eq. (46), one obtains the generalized GP equation [35] 



d 

ih— $(r,t) 



2 



|-V 2 + Kxt(r) + gn c (r, t) + 2gn{r, t) - iHR(r, t) 
2m 



(54) 



where 



27r\ /" dki c/k2 <ik 3 



(2tt) 3 (2tt) 3 (2tt) 3 
x 5(cj c + u;i - cj 2 - u; 3 )5(k c + ki - k 2 - k 3 ) 



x j[l + /(k l7 r, *)] [1 + /(k 2 , r, t)]/(k 3 , r, t) 

- /(ki, r, t)/(k 2 , r, *) [1 + /(k 3 , r, t)]}, (55) 

with p c = ftk c = ftV0, e c = /ic<j c = —hd t 9, and e(k;) = focuj being the condensate momen- 
tum, condensate energy and thermal cloud energy, respectively. We note that Eq. (54) is not 
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a closed equation because of the noncondensate mean field n(r, i) = (ip*(r,t)ip(r,t)) and the 
non-equilibrium distribution function /(k, r, t) in R(r, t). One should solve coupled equations 
that consist of the generalized GP equation and some kind of a kinetic equation, which is de- 
rived from the non-equilibrium Dyson equation (27), for the non-equilibrium distribution function 
/(k,r,*).[35] 

One can discuss collective modes of the condensate in the presence of the noncondensate atoms 
using only the generalized GP equation (54) as long as noncondensate atoms are in static equilib- 
rium, [55-57] or can be treated within the linear response theory. In the former case, the non- 
equilibrium distribution function in the dissipative term R(r, t) is replaced with the static equi- 
librium Bose distribution function. On the other hand, in the latter case, the dissipative term is 
neglected and only the mean field is considered by the linear response theory. In the previous 
paper, [49] we studied the microscopic mechanism of the Landau instability in a one-dimensional 
optical lattice using the generalized GP equation and the associated Bogoliubov-de Genne equa- 
tions derived from the Bose-Hubbard Hamiltonian with the static equilibrium approximation for 
the noncondensate atoms. 



B. Generalized GP Hydrodynamic Equations for the Condensate 

In order to study long-wavelength excitations, it is more convenient to use the hydrodynamic 
formulation. From the generalized GP action, Sggp^? G], which obtained in the previous sub- 
section, one can derive an action in terms of the density and the phase of the condensate order 
parameter by using the following variable transformation: 

$(r,t) = v /^M)e te ^, (56) 

where n c (r, t) and #(r, t) are the density and the phase of the condensate, respectively. The gen- 
eralized GP hydrodynamic action is given as follows: 



<Sggp [n c ,9,G] 
dr f dt 

I 



2 

h 2 



dt 



dt 



dr Jdt |-_^ c (r,t)V 2 v /n c (r,t) - ^n c (r,t) [V0(r,t)]' 
dvjdtn c (v,t) V CKt (r,t) + |n c (r,t) + 2gh(r,t) 
J drdv' J dtdt! yjn c (r, t)^n c {r', f) e -*[%.*)-% '•*')] j p( r> t ; r', if), (57) 
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where F is defined by Eq. (37). By minimizing the action (57) with respect to the density and 
phase and by performing the gradient expansion as performed in the previous section, one ob- 
tains the generalized GP hydrodynamic equations, which are equivalent to the generalized GP 
equation (54) [35, 55] 

dn c (r,t) 



ni- 



di 

d v c (r,t) 
dt 



V- [n c (r,t)v c (r,t)\ = -T 12 (r,t), 
m 



V 



0. 



(58) 
(59) 



where ri 2 (r, t) = 2n c (r, t)R(r, t) and the condensate chemical potential /i c (r, t) is given by 



h 2 V 2 Jn c (r,t) 

H c (r, t) = -- f + V ex t(r, t) + gn c (r, t) + 2gn(r, t). 



(60) 



2™ y/n c (T,t) 

The condensate velocity is defined by v = TiV6/m. The hydrodynamic equations (58) and (59) 
in the Thomas-Fermi approximation were used to discuss the damping of condensate collective 
oscillations in the harmonic trap potential at finite temperatures in Ref. 55. 



IV. COARSE-GRAINED FINITE-TEMPERATURE THEORY IN OPTICAL LATTICES 

Several authors have derived the hydrodynamic equations for the Bose condensate at zero tem- 
perature including the effect of a one-dimensional periodic lattice potential by focusing on the 
dynamics with length scale larger than a lattice spacing. [23, 24, 58] Using the hydrodynamic 
equations, Kramer et.al [24, 58] calculated frequencies of condensate collective oscillations. They 
have found that the frequency is renormalized through the effective mass due to the lattice poten- 
tial. [24, 58] Their results are found to be in good agreement with the experimental data obtained 
in Ref. 59. 

In addition to the experiments close to T = 0, such as in Refs. 5, 7, 8, 59, there have been very 
interesting experiments on the Bose condensate in optical lattices in the presence of the thermal 
cloud, such as damping in collective oscillation and the breakdown of superfluidity. [9, 20] Thus, it 
is tempting to derive finite-temperature hydrodynamic equations in the presence of the periodic lat- 
tice potential. In this section, we give a derivation of finite-temperature hydrodynamic equations, 
which are generalization of the hydrodynamic equations derived by Kamer et. al [24, 58] for the 
condensate, including the effects of the one-dimensional optical lattice at finite temperatures 
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A. Coarse-grained action 



In this subsection, we derive the coarse-grained action for Bose gases in optical lattice. Because 
the difficulty comes from the ^-direction trap potential, we first consider only a lattice potential in 
the ^-direction with a confining trap potential in the r_|_ directions. After that, the trap potential 
in the ^-direction can be included by the local density approximation. For this purpose, it is 
convenient to start with the generalized GP action r^[$, G] = Sqqp [*&, G], defined in Eq. (24), 
rather than the action defined in Eq. (57). Splitting this action into the three parts, one obtains 



S, 



GGP 



G] = S 1 + S 2 + S 3 



(61) 



where 



Si 



rd/2 
-d/2 

+ 2gn(r,t) 



S 2 



dz J dr ± Jjt<5>*{r,t) 
dz J dr ± J dt$*(r,t) 



^ d T r i \ q . ^ , . l9 

■^ + ^W + fl*('.*)l a 



9 h ^2 

-ih- 

dt 2m x 



Vhofrj 



$(r,t), 



S 3 = - j dzdz j dv ± dv' L j dtdt'<$>*(r,t)F(r,t)r\t')$(v',t'). 



(62) 

(63) 
(64) 



Here rj_ = (x,y) and Vj_ = d 2 / dx 2 + d 2 / dy 2 . V opt (z) is the optical lattice potential, while 14 Q (r) 
is the harmonic potential in Eq. (3). The function F(r, t; r', t') involved in the nonlocal part S 3 is 
defined by Eq. (37). The advantage of working with the action defined by Eqs.(62)-(64) is that the 
terms involve integrations over position and thus one can introduce coarse-grained approximation 
in a well-defined manner. We note that, in Eq. (62), the domain of the integration in the ^-direction 
is —d/2 < z < d/2 since we consider only the lattice potential in the ^-direction. 
For the variational parameter $(r, t), we use the following ansatz 



$(r,t) = $ fcc ( 2 )0 c (r ± ,t), 



(65) 



where &k c ( z ) = e u k c ( z ) * s me Bloch function with a quasi-momentum k c in a lowest band. 
The condensate Bloch amplitude Uk c (z) has the periodicity of the lattice potential. This ansatz is 
= and Vh (z) = 0, i.e., uniform in the ^-direction. Similarly, the field operator 



exact for V c 



opt \ 



for the noncondensate is assumed to be expanded by the Bloch states: 

n k 



(66) 
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where y3 n ,fc(r_i_) is a destruction operator of the noncondensate in the r^-direction. The Bloch 
function <j> n ,k{z) — e lkz u n ,k(z) describing the thermal cloud atoms with a quasi-momentum in a 
n-band, which satisfies the following Schrodinger equation 

H{z)4> n , k (z) = 4°l0n,fc(>), 



where 



H(z) 



h 2 d 2 



+ V opt (z). 



2m dz 2 

The Bloch amplitudes of both condensate and thermal atoms are assumed orthonormal, 



fd/2 

/ dz u* kc (z)u K (z) = S kciK , 

J — d/2 



rd/2 



dz u* k {z)u nl>k '{z) = 5 n>n '8 kik '. 

J -d/2 

Substituting Eqs. (65) and (66) into Eq. (62), one obtains 

r d / 2 r r r h 2 d 2 

Si = - I dz dr ± dt$* kc {z)cp* c {r L ,t) 



-d/2 

+ §l^)| 2 |^(r ± ,t)| 2 



2m~d? + Vopt[Z 



+ 2 9^2Y1 C,fc(^)0n',fc'( 2; )(^l,fc( r ± ; t)fi n/)k/ {r ± , t)) 

n,n' k/k' 

-J dr± J^dt |v9 c (r_L,t)| 2 



®k c (z)<-Pc(r±,t) 



rd/2 

x / , dz$* kc (z) 

J-d/2 



+ 2^EE K,k(z)<t>n>,k>{z)((pl k {r ± , t)<p n > jk >(r±, t)) 

n,n' k,k' 

Using the Bloch amplitude u kc (z) and u k (z), Eq. (71) becomes 

St = - J dr ± J dt \ip c (r ± ,t)\ 2 

rd/2 r *2 



x / , dzul (z) 

J-d/2 

9 



+ ||wfc c (^)| 2 |^(r±,t)| 2 
+ 2 9 H K,k( z )^n',k'(z)(0l tk (r L , z, t)0 n ,, k ,{r ± , z, t)) 

n,n' k,k' 

where <p k (r±, z, t) = e tkz (p k (r±, t). Similarly, Eq. (63) becomes 

r a t2 



S, = 



dr± J dt y?*(r_|_, t) 



(67) 

(68) 

(69) 
(70) 



(71) 



(72) 



(73) 
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where we have used the normalization conditions, Eqs. (69) and (70). 

Finally we consider S3. In the system only with the harmonic potential, the Wigner transforma- 
tion was performed by expanding the Green's function by a plane wave, as reviews in Sec. III. On 
the other hand, in the system only with the one-dimensional periodic lattice potential, the Green's 
function should be expanded by the Bloch function, which has the periodicity of the lattice poten- 
tial. The Green's function expanded by the Bloch function is given by 

iG{v, t; r', f) = E E hk(z)fa,k>( z ') < T c[vMr±, t)0l jk ,(v' ± , t')}) 

n,n' k,k' 

= EE^,^)^,.(^e ! ^" lW (Tc^, i (ri,i)^ ) ,(r' 1) t')]). (74) 

n,n' k,k' 

The Bloch amplitudes u n>k (z) and u* n ^ k {z) satisfy Eq. (67). It will be convenient to define the 
following Green's functions: 

z^/(r ± ,t;r' ± ,t') = (T c [0 n , fc (r ± , t)^, fc ,(r' ± , *')]>, (75) 
iG n k $(r ± ,z,t;r' ± ,z',t') = e ik *e- ik '*' (T c [0 n , k (r ± , t)^^, if))). (76) 

In terms of g2' k > or G k ' k , , the Green's function G(r, t; r', t') can be written as 

G(v,t;v',t') = ( 77 ) 

n,n' k,k' 

= EE <k(z)K'A z ') Gn kS'( r ^ z > *; r '±> z '> 0- ( 7g ) 

n,n' A;, ft' 

These Green's functions are useful when we apply the coarse-graining procedure. With use of 
Eqs. (65), (75), and (76) in (64), one obtains 

5 3 = -yEEEEEE 

ni,nj ri2,n' 2 713, rig fci,fc^ fc2,fc 2 k3,k' 3 
rd/2 „ 

x / d^ ni)fel (z)0 n2)fe2 (^)^, fe ,(z)$^ (z) 

/•d/2 _ 
J-d/2 11 2 2 

x /" dr ± dr' ± J dtdt' ip* c (r ± ,t)ip c (r' ± ,t') 

x (r-u *5 r'x, 0^(r±, *5 ^, i')<Of ft, f ; r ± , t). (79) 

Now we include a confining trap potential in the ^-direction. As in the zero-temperature case, 
we are only interested in the dynamics on a length scale much longer than the lattice constant d, 
and thus the local density approximation can be used to treat the trap potential. In addition, we 
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assume that one can use the Bloch function for describing the condensate on this length scale. 
We assume that the axial size of the condensate, R z , in the confining trap potential is much larger 
than the lattice constant d. The difficulty comes from the rapid variation of the length scale by 
the lattice potential. In order to treat this length scales appropriately, we first regard the terms 
e lkcZ if c (r±,t) and e lkz (p njk (^±,t) in Eqs. (65) and (66) as almost constant in the length scale of 
the lattice constant —d/2 < z < d/2, although, which still depends on z. Then the variation due 
to the length scale of the confining trap potential occurs in the length scale much longer than the 
lattice constant. The condensate order parameter and the Green's function for the noncondensate 
are assumed to vary as the site index I: 

$ c (r,t) - u kc (zy k ^ z v c (v ± ,l,t) 

= u kc (z)$ c (r ± ,l,t), (80) 

G(r,t-y,t>) -> EE^(^<^> lt ' (riil)2 ^ (r - iK 

n,n' k,k' 

x g k ,k> ( r ±,M; r±,i,t) 

= Y.l^UnA z )K>,k>( z ')Gk t t l {r ± J,t-,r' ± ,l',t'). (81) 

n,n' k,k' 

We can also define the site-represented amplitude and phase of the condensate order parameter: 

$c(r_L, /, t) = ^n c (r ± ,My s ^), (82) 

where 

n c (r ± ,l,t) = |$ c (r ± ,/,t)| 2 , (83) 
S(r ± ,l,t) = k c (r ± ,l t t)z + 0(T ± ,l,t). (84) 

In Eq. (84), ??(rj_, /, t) is the phase of (p c (r±, /, t). By using Eq. (84), a site-represented condensate 
velocity can be defined by 

h d 

v z (r±,l,t) = ——S(r ± ,l,t) 

m oz 

= -k c (r x ,l,t), (85) 
m 

v x Jr ± ,l,t) = -V x S(r ± ,l,t). (86) 
m 

The above site-represented quantities can be understood by identifying these quantities as averaged 
one: 

r ld+d/2 
-d/2 



I Ad+d/2 

n c (r±,l,t) = - / dzn c (r,t), (87) 
d Jid-d/2 
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X Ad+d/2 % Q 

v c (r ± , I t) = - / dz — — 0(r, t), 

d Jld-d/2 m oz 



(88) 



1 fld+d/2 ^ % Q 
-d/2 

where n c (r, t) = |$(r, t)\ 2 and 9(r, t) are the condensate density and phase, respectively, which 
are the solution of the generalized GP equation (54). One can also introduce the averaged Green's 
function 



I Ad+d/2 



dz 



Vd+d/2 



dz' G^,(r,t;r',0 



(89) 



d 2 Jld-d/2 Jl'd-d/2 

where G k ' k , (r, t; r', t) is the solution of some kind of kinetic equation for the noncondensate atoms, 
for instance, the Kadanoff-Baym equations [54]. 

Substituting Eqs. (80) and (81) into Eqs. (72), (73), and (79), one obtains 

rd/2 



Si = -£/ dr ± J c dt\<S> c (r ± J,t)\ 2 J_^dzul(z) 



x 



d 



2^ ( Tz + ik? ) + V opt (z) + | \u kc (z)\ 2 \<p e (r ± , Z, t)f 



2 # E E <,k( z ^n>,k> (z)iGl%, < (r ± , Z, t; r ± , I, t) 

n,n' k,k' 



(90) 



S2 = —^2 J dr± J dt 



x*;(rj.,(,i) 



$ c (r_L,Z,t), 



(91) 



S, 



EE E E E E E 



nijTl^ 712, ^2 n 3! n 3 ^li^i ^21^2 ^31^3 

x y dr ± dr' ± Jjtdt' $*(r ± , Z, i)$ c (r' ± , Z', f) 
G k l'^(r-±, I, t; r'j_, I', t')G £ k ?(r ± , Z, t; r' ± , I', t')G k *' k ?(r' ± , I', t'; r ± , Z, t) 



x 



x 



d/2 

-d/2 
d/2 



dz U 



niM\ Z ) U n2M\ Z ) U n'k'\ Z ) U k c \ Z ) 



x / , dz' u* n , k , (z')u* n/ k , (z')u n3M {z')u kc {z'), 

J -d/2 11 2 2 



(92) 



where Vh (r±, Z) is a confining harmonic trap potential which labeled by Z in the z-direction. 

The averaged quantities in Eqs. (90), (91), and (92) are assumed to be smooth functions of r_|_ 
and to vary slowly with the lattice site index Z. Thus, the coarse-grained macroscopic densities and 
velocity can be obtained by replacing the discrete index Z with the continuous variable z = Id. This 
retains only the information on length scale much larger than the lattice spacing d by focusing on 
the "macroscopic" dynamics. Thus the "microscopic" information shorter than the lattice constant 
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is averaged out and only enters the modified local condensate energy given below. We define the 
coarse-grained quantities as follows: 



$ c (r±,Z,t) 
n c (r ± ,l,t) 
v c (r±,l,t) 



$c{r±,z = ld,t), 
n c (r±, z = ld,t), 
v c (r_L, z = ld,t), 

G k ' k , ( r -L' z = r J-' z ' = I'd, t')- 



The coarse-grained phase of the condensate is related with u c (r, t) by the following equation: 

VS(r,t) = jv c {r,t). 

Note that this equation gives the definition of the coarse-grained phase of the condensate. 
With the coarse-grained quantities, we can obtain the coarse-grained action 



(93) 
(94) 
(95) 
(96) 

(97) 



S, 



CG 



S± + S2 + S3 



where 



Si 



J dr J dt $ c (r,t) e opt [k c ,r,t], 

d h 2 



S 2 = - dr dt$*{r,t) 



S 3 = - J drdr' J dtdt'<f>* c (r,t)F(r,t;r',t')<S> c (r',t'). 



Non local term F(r, t; r', t') in S3 is given by 



F{r,t;r',t' 



yEEEEEE 

rai,rij n2,n' 2 713, rig k2,k' 2 k2,k' 2 

x (r, t; r>, t')G^ (r, t; r', t')G^ (r', f; r, t) 

[ d l 2 

x / dz u nukl (z)u n2>k2 (z)u* n , k , (z)u* kc (z) 

J-d/2 3 3 

[d/2 

X / , dz u* n , k > (z)u* n , Az)u n3M (z')u kc {z'). 

J-d/2 11 2 2 



In Eq. (99), we have defined 



^opt [k c , r, t] 



d/2 
-d/2 



dz u* kc (z) 



2m \ dz 
1 1 



+ ik c 



V opt (z) 



u kc {z) 



g ra/z 

+ TT^cM) / ; \u kc {Z 
2 J-d/2 

/ari 
dz \u kc (z)\ 2 <, 

n,n' ' 



k (z)u n ,, k ,(z) 



(98) 

(99) 
(100) 
(101) 



(102) 



(103) 



22 



B. Coarse-grained generalized GP equation 



In this subsection, we derive the coarse-grained GP equation from Eq. (98). By taking the 
derivative with respect to <£>*, one obtains 

ih^$ e (r,t) = |-^Vi + e opt (A; c ,r,t) + y ho (r)|<l c (r,t) 

+ J dr' j dt'F(r,t;r',t')$ c (r',t'). (104) 

This equation involves the effects of the lattice potential by e opt (k c , r, t) and of the coupling to the 
thermal cloud by the non-local function F(r, t; r', t'). The crucial point here is that the difficulty 
which comes from the rapid variation due to the lattice potential disappeared and the effects of the 
lattice potential is effectively included. 

In order to further simplify the correlation function F, we assume that the coarse-grained 
macroscopic variables vary slowly in space and time compared to the spatial and temporal scale 
of a collision event. We can then approximate the condensate order parameter at (r', t') close to 
(r, t) by a Taylor expansion 



$(r',t') = Jn c {r>,t')e iS{r '^ 



n c (r,t)e i[5(r ' i)+at5{r '* )( *'~' )+v,s(r ' <) ' (r '~ r)1 
= <5(r, £) e ~ i M r ' t )(*'~')~ kc ( r '*H r ' -r )]. (105) 

The condensate frequency and wavevector are defined by uj c {r,t) = —d t S(r,t) and k c (r, t) = 
VS(r,t), respectively. For the Green's function, we follow the same procedure as in Sec. Ill A 
by using the Wigner transform. We note, however, that the Green's function Gfc zl ,fc z2 (ri, t\\ r 2 , t 2 ) 
in this section involves the band index n and quasi-momentum k z . We first expand the Green's 
function by the plane-wave and rewrite by the center-of-mass, R = (ri + r 2 )/2, and relative 
coordinates, r = r x — r 2 



r n,n'< , , * \ . . Jkziz-L -ik z2 z 2 f dk ±A f ^ k -L,2 



(2n) 2 J (2tt) 2 
x e ik±i r±1 e- ikx2 r±2 G n ' n ' < (k 1 , k 2 , h, t 2 ) 



J n,n 



(2vr) 2 J (2tt) 2 

xe* 2 V kl - Rl G n < (k ) K,M'), (106) 

where we have introduced K = (ki + k 2 )/2 and k = ki — k 2 and neglect the multi-band effects in 
the Green's function. We have used a notation: G<(k, K, t, t') = G n - n< (k, K, t, t'). The Wigner 
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transformed Green's function is defined by 

duo 



G<(K>;R,T) = J £e^^G<(K,K,t,t') 



duo 



-iui(t~t') 



dk 



e lk -*G<(k,K,t,t'), 



(107) 



2ir~ J (2tt) 3 

where the integral for the quasi-momentum k z comes from the one in Eq. (102). We then use the 



quasi-particle approximation, changing the notation K — ► k 

iG<(k, uo; R, T) = 2tt5(u - e n (k; R, T)/h)f n (k, R, T), 
iG>(k, uo; R, T) = 2<k8{u - e n (k; R, T)/h) [1 + f n (k, R, T)], 



(108) 
(109) 



where the noncondensate energy e is derived from the equation of motion for the retarded Green's 
function and is given by 

h 2 k 2 t ^ . 



e n (k;R,T) = —-^ + e^ kz + 2gn c (R, T) / dz \u kc {z)\ 2 \u n ^{z 



d/2 
d/2 



Um,qX Z )\ 2 \^n,kX Z 



(110) 



where ef ) kz is a solution of Eq. (67). With use of the Wigner transformed Green's function (107) 
and the quasi-particle approximation (108), Eq. (103) becomes 



Uu (Z 



e opt (k c ,R,T) = J ^dz 

+ -n c (R,T) / K c (2:)| 4 

/• dk r d / 2 
+ 2 #E / 77T^/«( k ' R ' T ) / dz \u n ,k z 



»l a K(*)l 2 - 



(111) 



With these approximations, we can derive the generalized GP equation in an optical lattice 
potential, replacing (R, T) with (r, t) 

i%— $ c (r, t) = {- — Vi + e opt {k c , r, t) + V ho (r) + iR(r, t) }<l c (r, *). (112) 



where the dissipative term R(r, t) is given by 

•2tt^ 4 



R(r,t) 



x 



E E 

mGN ni,n2,n,3 
d/2 

dzul 

-d/2 



dk\ dk-2 dk% 
(2tt) 3 (2tt) 3 (2tt) 3 
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x 5(k ±c + k ± i - k ±2 - kj. 3 ) 

x 5(k zc + k z i - k z2 - k z3 - 2mq B ) 

x | [1 + f ni (k a , r, t)} [1 + f n2 (k 2 , r, t)} f n , (kg, r, t) 

- f ni (k l5 r, t)/„ 2 (k 2 , r, t) [1 + /„ 3 (kg, r, t)] }. (113) 

Here p c = hk c is the condensate momentum and qs = 7r/d in the umklapp term Imqs (where 
m is an integer) in the 5-function is the Bragg wave-number . Compared Eq. (1 13) with Eq. (55), 
we see two new features associated with the lattice potential. First, the bare coupling constant g is 
renormalized by the Bloch functions. Secondly, the momentum conservation for the z-direction is 
modified to the looser condition: Tik Z)C + %k Zj \ — hk z>2 — Tik z ^ = 2mhqB- This reflects the break- 
down of the translational symmetry due to the periodic optical lattice potential. In the formalism 
developed here, the collision dynamics information of a length scale shorter than the lattice spacing 
d is effectively included through these two new features. The local condensate energy, e c = Huj c , 
is defined by 

e c (r,t) = - V± V^' j) + (fc ^ + Ho(r) + ™^l ± , (114) 
2m ^n c (r,t) 2 

where /i op t(^c, r, t) = d(n c e opt ) /dn c . This expression for the condensate energy is given below 
after discussing the quantum hydrodynamic formulation for the condensate by the coarse-grained 
quantities. The noncondensate energy, q = Tiuji, in Eq. (113) is defined by e« = e ni (kj, r, t). 



C. Coarse-grained generalized GP hydrodynamic equations 

In order to derive the coarse-grained equations in terms of hydrodynamic variables, we should 
start from the coarse-grained action because e opt in the coarse-grained generalized GP equa- 
tion (112) depends on the condensate velocity through the wavevector k c , preventing usual vari- 
able transformation from <5 C and $* to n c and v c in the coarse-grained generalized GP equation 
Eq. (112). 

Combining the results given in Eqs. (99), (100), and (101) with Eq. (103), and using the hydro- 
dynamic variables, Eqs. (94), (97), we can write down the effective action for the condensate in 
optical lattices in terms of the coarse-grained macroscopic variables 



Sefl: = I dr I dt 



c 



ihdn c (r,t) dS(r,t) 
- hn c (r,t) 



dt cx ' ' dt 
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(115) 



The quantum hydrodynamic equations for the condensate are found by minimizing the effective 
action (115) with respect to the density and phase, leading to two coupled equations: 

h 2 Viv^eM) 
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m — h V 
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dk. 
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n c (r,t) 

Ar,t)e- iS{r ' l) I dr' I dt'F(r,t;r',t')Jn c (r',t')e iS( - T '^ 



where we defined the condensate chemical potential 

d (n c E pt(k c , n c , n nc )) 



(117) 



(118) 



The condensate equation of motion given by Eqs. (116) and (117) describe the long-wavelength 
dynamics of the condensate in the presence of the periodic lattice potential, including the crucial 
coupling to the thermal cloud described by the function F(r, t; r', t'). All information associated 
with the length scale shorter than the lattice spacing d is buried in the Bloch functions. 

Similarly to the coarse-grained generalized GP equation, one can simplify the above equations 
by local approximations and obtain the following equations: 

h 2 VWn c (r,t) 
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2m 
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-f(r,0, 

(120) 

where the source term in Eq. (120) is given by 

f(r,t) = 2n e (r,t)R(r,t), (121) 

where -R(r, t) is defined by Eq. (113). 

At zero temperature, the GP hydrodynamic equations for the lattice potential have been derived 
by the tight binding approximation in Ref. 24, where the equation of state for /x opt is assumed to 
have the same structure of the uniform system. Ref. 23, 26 have obtained the GP hydrodynamic 
equations for a general /i opt and e opt . Compared with the GP hydrodynamic equations derived in 
the above works, Eqs. (119) and (120) are finite-temperature generalization, which includes the 
interaction between the condensate and noncondensate atoms. At T = 0, of course, n in /x opt and 
e opt and f vanishes, and thus Eqs. (119) and (120) reduce to the GP hydrodynamic equations for 
the lattice potential. 



1. Low condensate velocity approximation 

In this subsection, we restrict ourselves to a condensate moving with a small superfluid velocity. 
In that case, it is very useful to introduce the effective masses for the long wavelength in the lowest 
Bloch energy and chemical potential band and rewrite Eqs. (119) and (120) by using these effective 
mass. The effective mass and chemical potential effective mass are defined by the following 
equations [23, 24, 26]: 

1 . d 2 e opt (k c 
ml ~ h 2 dk 2 c 

1 _ 0W(*c, (123) 



(122) 

fc c =o 



m* i0 tfdkl 

In the usual Bloch theory of metals, we have only one effective mass, namely, rrtg = Tt* . For 
the Bose condensate, however, due to the interaction term the two different effective masses enter 
the theory. This is originally comes from the two energies, the condensate energy band e opt and 
chemical potential band ji opt . In terms of these effective masses, (122) and (123), the coarse- 
grained generalized GP hydrodynamic equations, given by (1 19) and (120), become 
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rrin 



= -r(r,t), (125) 

where /i op t(^e, n nc ) = [i op t(k c = 0, n c , n nc ). In the low condensate velocity approximation, the 
local condensate energy (1 14) is given by 



e c (r,f) 
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(126) 



The local condensate chemical potential is given by 

// c (r,t) = fXopt(n c ,n nc ) + V ho (r). (127) 

When we set T — > and /i op t(^c, n n c) — ► ^ op t(n c , n nc = 0), appropriate when the thermal cloud is 
absent, Eqs. (124) and (125) reduce to those obtained in Ref. 24. 



V. INSTABILITY OF THE SUPERFLUID 



The Landau instability of the Bose condensate has been studied by using the original GP equa- 
tion. [21-26] Within the GP equation, one can derive the stability phase diagram for the condensate 
from the negative excitation energy condition. However, this approach does not give any insight 
into the microscopic origin of the Landau instability. Moreover, one has to construct a micro- 
scopic theory for the landau instability because the original Landau argument cannot be applied 
to the lattice system where the momentum conservation in not satisfied, which Landau originally 
relied on. 

In the present section, we use the finite-temperature theory developed in the previous sections 
to exhibit a specific microscopic origin of the Landau instability of superfluidity in a Bose con- 
densate. 



A. Relation between the damping of collective modes and the instability of the condensate 

We shall show that f defined by Eq. (121) can change sign and as a result leads to a Landau 
instability. [49-51] This is illustrated generally as follows. The amplitude of the collective mode 
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S$ Pct k decays as 8$ Pc> k oc e~ Fpc - kt , where r Pc> k is a damping rate. This relation indicates that 
the superfiuid state is stable as long as T Pc k is positive for any collective mode momentum k. In 
fact, for a condensate at rest (p c = 0) condensate, one can show that the damping rate is always 
positive, [49, 50] and thus the collective mode decays exponentially in time. However, as shown in 
the following subsections, f Pc can become negative in the case of a moving condensate (finite p c ). 
A negative value of T Pc k indicates an instability. The exponential growth in time of the amplitude 
of collective modes due to collisional coupling to the thermal cloud destabilize the condensate 
state, leading to the breakdown of superfluidity. Recent experiment [9] appear to support this 
scenario. 

In the collisionless regime of interest to which we restrict ourselves, there are two important 
damping processes, namely, collisional damping and Landau damping. The former occurs due to 
the collisional exchange of atoms between the condensate and noncondensate, [55, 56] while the 
latter occurs due to the dynamical coupling between the condensate oscillation and the thermal 
excitations. [40-42, 60-62] In the present paper, we will only consider collisions between the 
condensate and noncondensate atoms, and also ignore the harmonic trap potential (Vh (r) = 0). 

B. Instability due to the Collisional Damping Process 

The collisional damping process considered in the present paper was investigated by Williams 
and Griffin [55, 56] and Duine and Stoof [57] for collective modes of the condensate in the har- 
monic trap potential. The collisional damping arises due to the lack of diffusive equilibrium be- 
tween the condensate and noncondensate, namely, equilibration process due to the collisional ex- 
change between the condensate and the noncondensate. The authors of Refs. 55-57 calculated the 
collisional damping rate by solving the dynamical equation of motion for the condensate, while 
the noncondensate is treated as being static. To simplify the notation, we omit the bars in this 
section, however, we note that n c , v c , and T always refer to coarse-grained values. 

We approximate the non-equilibrium distribution function in Eq. (1 13) by the static Bose dis- 
tribution function for thermal equilibrium [55, 56, 63] 



where JIq is the chemical potential of the noncondensate atoms and the energy of the noncon- 
densate atoms is given by the Hartree-Fock approximation appropriate to one-dimensioal lattice 




(128) 
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potential along the z-axis; e„(k) = hu n (k), which is defined by Eq. (110). Our use of the static 
thermal cloud approximation implicitly assumes that the thermal excitations reach equilibrium 
with a relaxation time much shorter than the period of the condensate collective modes. This as- 
sumption may be justified by the experiment [20], where the thermal cloud in the lattice potential 
reaches its equilibrium state very rapidly. 
Using the identity for the Bose distribution 

A(l + / 2 )(1 + / 3 ) = (1 + /i)/ 2 / 3 e^-^-^-«), (129) 

the source term T(r, t) in Eq. (121) reduce to 
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We will derive linearized equations of Eqs. (119) and (120) with n c (r,t) = m? c + 6n(z) and 
v c (r, t) = [v ® c +5v c (z)]z, where n° c and v° zc are static values of the condensate density and velocity, 
respectively. We need to evaluate V to first order in the deviations away from equilibrium. The 
condensate energy (126) is expanded as 



^ + m wk- c 5Vzc + ^ 6nc 

/zo + m 5v zc H — — — ora c . (131) 



where /i° pt = fi op t(k c0 , n c0 , n nc ,o), and n c0 and v c0 = ft/c c o are the condensate density and velocity 
in equilibrium, respectively. Here we have used the relation fjP t = fl . Note that the fluctuation 
of the noncondensate is neglected because we use the static thermal cloud approximation. 
Using Eq. (131), one can also expand e^ ec ~^ in Eq. (121) to give 

exp[/3(e c - /2o)] 
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exp < p 
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The resulting linearized form of the dissipation term T(r, t) reduces to 

5T(z,t) 



^—Sn c (z, t) + m-^Sv^z, t) 



(133) 



dn c ' hdk c 

where the relaxation time r arising from the collisions between the condensate and noncondensate 
atoms is defined by 

dki dk 2 dh 3 
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where the superscript "0" of uj® and oj® (where i = 1,2,3) indicates the quantities take its static 
value. We use Eq. (133) in the linearized version of the generalized GP hydrodynamic equations, 
leading to 
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(135) 
(136) 



where we have defined the effective mass and the chemical potential group velocity at the arbitrary 
value of the condensate velocity k° c as follows [24, 25]: 

1 d 2 e r 
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The effective mass m* is a finite- k c generalization of for the long wavelength defined by 
Eq. (122). We note that the chemical potential group velocity v% is deferent from the usual group 
velocity because v% is derived from the chemical potential band, while the usual group velocity is 
derived from the energy band e opt . 

Before solving the coupled equations (135) and (136) for the condensate fluctuations, it is 
useful to derive the Stringari-type equation for the condensate fluctuation [64], which gives the 
frequency of the condensate collective modes, in order to show the significance of the collisions 
between the condensate and noncondensate atoms buried in r. One can show that the linearized 
equations (135) and (136) reduce to the finite-temperature Stringari equation in the presence of the 
lattice potential when we set v% = 0: 



d 2 Sn c (z,t) n° c d 2 



opt 



dn c 



Sn c (z,t) 



r' at ' K ' 



(140) 



dt 2 m* dz 2 

where 

i _ /K<>t 

t' r dn c 

The finite-temperature Stringari equations was first derived by Williams and Griffin [55] for a 
harmonic potential. The collision time r' describes collisions between the condensate and non- 
condensate atoms when the condensate is perturbed away from equilibrium. Eq. (139) clearly 
shows that the new term associated with the collision time r' in Eq. (139) causes damping of the 
condensate fluctuations. This damping is due to the lack of collisional detailed-balance between 
the condensate and the static thermal cloud pinned by the lattice potential. The collisional damping 
is important damping process in addition to the Landau damping in the collisionless regime. 

To solve the coupled equations (135) and (136), we assume a plane-wave solution ~ 
exp[i(q z z — ut)] for both 6n c and Sv c , and then Eqs. (135) and (136) give 

in 1 / 9/jP \ 

uj8n c -n° c q z 5v c - v^q z 5n c = i- mv^5v c + ° pt 5n c , (141) 

m* t V on c ) 

1 dfjP 
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m on c 

From Eq. (142), one finds 

5v c = —^^qz ^—rr- 5n c . (143) 

m on c to — Vcq z 

Substituting this into Eq. (141) and eliminating Sv c , one obtains 

-^^r-^ = V w - (144) 
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In the absence of the collisions, i.e., taking the limit 1/V — > 0, we find the collective mode 
frequency f2 given by [23-26, 58] 

u = v»q z ± c*q z = a (145) 

Here 



c' = J4^, (146) 
V m* on c 

is a Bogoliubov-type sound velocity modified due to the presence of the lattice potential, while 
v% is defined by Eq. (138). In Eq. (145), the opposite sign "±" correspond to a sound wave 
propagating in the same and in the opposite direction, respectively. Kramer et al. [58] gave the 
sound velocity for an optical lattice by using tight-binding model with a specific approximation 
for the chemical potential. The general expression of the sound velocity was given by the GP 
hydrodynamic analysis by Machholm et al. [23] and by Kramer et al. [24]. Taylor and Zaremba 
used the Bogoliubov equation by a systematic expansion in powers of the phonon wave vector [26] . 
Our result of the sound velocity is natural extension of above works at zero temperature to finite 
temperature. For a translationally invariant system and at zero temperature, /i op t — > gn° c and 
m* — > m. Therefore, the sound velocity is given by the usual result c = J 

Without the collision term, the effect of the thermal cloud enter into the collective mode fre- 
quency by the mean-field interaction in the chemical potential. In this case, the collective mode 
does not damp. Now we study the effects of the collisions between the condensate and noncon- 
densate atoms represented by the collision term t', which is second order effect of the coupling 
constant. For this purpose, we include the effects of the collisions represented by the relaxation 
rate 1/V, giving the dispersion relation as w = - iT c . To first order in 1/V, one obtains the 
collisional damping rate to be 

r c = — (l±-5- . (147) 



2V V c* 

This is the key relation to consider the instability of the condensate. Recall the argument in the 
first part of this section. As far as the damping rate T c is positive, the collective mode of the 
condensate is stabilized by this damping process. This corresponds to the case of lower sign in 
Eq. (147). This expression, however, shows that T c can be negative, indicating a growth instability 
when the direction of the condensate and sound velocity is opposite, and when 

v» > c*. (148) 
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This condition turns out to be same as the usual Landau criterion for the superfluidity in a uniform 
system, except that the condensate and the sound velocities are now modified due to the presence 
of the optical lattice potential. This type of instability has been discussed by imposing on the 
condition that the excitation energy becomes negative, i.e., hCl < 0. The region of this instability, 
so called Landau instability, for an optical lattice potential was first obtained by Wu and Niu [21, 
22]. The same argument was performed by several authors [23, 25, 26, 28]. In contrast to the 
previous works, the crucial point in the present work, however, is that we derive the Landau 
criterion by specifying the microscopic destabilization process, which is performed by calculating 
the damping rate of the condensate collective mode. In this sense, we give a explanation of the 
microscopic mechanism of the Landau instability. This kind of discussion for the stability can be 
also seen in Refs. 63, 65 for a trap potential and uniform system, respectively. 

VI. SUMMARY AND CONCLUSIONS 

In this paper, we have developed a coarse-grained finite-temperature theory for a Bose conden- 
sate in one-dimensional optical lattices, in addition to the confining harmonic trap potential. This 
theory consists of coarse-grained equations of motion for the condensate variables and noncon- 
densate Green's functions, which include the effect of a dissipative term due to collisions between 
the condensate and the thermal cloud, as well as the noncondensate mean-field. 

With use of the non-equilibrium field theory, the 2PI effective action for the Bose condensate 
on the Sch winger- Keldysh closed-time path has been obtained. Introducing an ansatz for the 
variational function in the effective action to perform a coarse-graining approximation, we have 
obtained a coarse-grained effective action, which includes the effects of the optical lattice potential 
effectively, in the presence of a thermal cloud of noncondensate atoms. We have also derived a 
coarse-grained action in terms of hydrodynamic variables of the condensate. Using the variational 
principle, we obtained coarse-grained equations of motion for the condensate variables, which can 
be used to describe the long wave-length dynamics on the length scale much longer than the lattice 
constant d. 

To illustrate our formalism, we used the generalized GP hydrodynamic equations to investigate 
the stability of superfluidity in the current-carrying condensate. Following recent work [49, 50], 
we calculated the damping rate of the collective oscillations. We have found that the collisional 
damping rates change sign when the condensate velocity exceeds the renormalized sound velocity, 
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leading to the Landau instability consistent with the Landau criterion. The results in this paper 
sheds light on the microscopic origin of the Landau instability. 

In the present paper, we concentrated on the effect of the optical lattice and ignored the trapping 
potential. One could use our formalism to analyze the experimental results by Florence group [5, 
20] on the damping of the condensate in dipole oscillations due to the thermal cloud. 

In order to describe the coupled non-equilibrium dynamics of both the condensate and non- 
condensate, one has to derive a kinetic equation for the noncondensate distribution function in the 
presence of an optical lattice. The generalized GP hydrodynamic equations for the condensate 
derived in the present paper and the kinetic equation for the noncondensate will be used as a sound 
basis for investigating finite-temperature behaviors of the Bose condensate in optical lattices. 

ACKNOWLEDGMENTS 

The authors thank Allan Griffin for useful suggestions and a critical reading of the manuscript. 
The authors also acknowledge valuable and helpful discussions with Ippei Danshita and Kiyohide 
Iigaya. S. K. is supported by JSPS (Japan Society for the Promotion of Science) Research Fellow- 
ship for Young Scientists. 



[1] I. Bloch, Nature Phys. 1, 23 (2005). 

[2] O. Morsch and M. K. Oberthaler, Rev. Mod. Phys. 78, 179 (2006). 

[3] M. Greiner, O. Mandel, T. Esslinger, T. W. Hansen, and I. Bloch, Nature 415, 39 (2002). 

[4] K. Xu, Y. Liu, J. R. Abo-Shaeer, T. Mukaiyama, J. K. Chin, D. E. Miller, W. Ketterle, K. M. Jones, 

and E. Tiesinga, Phys. Rev. A 72, 043604 (2005). 
[5] S. Burger, F. S. Cataliotti, C. Fort, F Minardi, and M. Inguscio, Phys. Rev. Lett. 86, 4447 (2001). 
[6] F S. Cataliotti, L. Fallani, F. Ferlaino, C. Fort, P. Maddaloni, and M. Inguscio, New. J. Phys. 5, 71 

(2003). 

[7] L. Fallani, L. De Sarlo, J. E. Lye, M. Modugno, R. Saers, C. Fort, and M. Inguscio, Phys. Rev. Lett. 
93, 140406 (2004). 

[8] C. D. Fertig, K. M. O'Hara, J. H. Huckans, S. L. Rolston, W. D. Phillips, and J. V. Porto, Phys. Rev. 
Lett. 94, 120403 (2005). 



35 



[9] L. De Sarlo, L. Fallani, J. E. Lye, M. Modugno, R. Saers, C. Fort, and M. Inguscio, Phys. Rev. A 72, 
013603 (2005). 

[10] M. Albiez, R. Gati, J. Foiling, S. Hunsmann, M. Cristiani, and M. K. Oberthaler, Phys. Rev. Lett. 95, 

010402 (2005). 
[11] D. Jaksch and P. Zoller, Ann. Phys. 315, 52 (2005). 
[12] E. Altman and A. Auerbach, Phys. Rev. Lett. 89, 250404 (2002). 

[13] A. M. Rey, B. L. Hu, E. Calzetta, A. Roura, and C. W. Clark, Phys. Rev. A 69, 033610 (2004). 
[14] A. M. Rey, B. L. Hu, E. Calzetta, and C. W. Clark, Phys. Rev. A 72, 023604 (2005). 
[15] E. Altman, A. Polkovnikov, E. Demler, B. I. Halperin, and M. D. Lukin, Phys. Rev. Lett. 95, 020402 
(2005). 

[16] C. Kollath, A. M. Lauchli, and E. Altman, Phys. Rev. Lett. 98, 180601 (2007). 
[17] K. Temme and T. Gasenzer, Phys. Rev. A 74, 053603 (2006). 
[18] L. D. Landau, Zh. Eksp. Theo. Fiz. 11, 592 (1941). 

[19] A. A. Abrikosov, L. P. Gorkov, and I. E. Dzyaloshinski, Methods of Quantum Field Theory in Statis- 
tical Physics (Dover, New York, 1975). 

[20] F. Ferlaino, P. Maddaloni, S. Burger, F. S. Cataliotti, C. Fort, M. Modugno, and M. Inguscio, Phys. 
Rev. A 66, 011604(R) (2002). 

[21] B. Wu and Q. Niu, Phys. Rev. A 64, 061603(R) (2001). 

[22] B. Wu and Q. Niu, New J. Phys. 5, 104 (2003). 

[23] M. Machholm, C. J. Pethick, and H. Smith, Phys. Rev. A 67, 053613 (2003). 

[24] M. Kramer, C. Menotti, L. Pitaevskii, and S. Stringari, Eur. Phys. J. D 27, 247 (2003). 

[25] C. Menotti, A. Smerzi, and A. Trombettoni, New J. Phys. 5, 1 12 (2003). 

[26] E. Taylor and E. Zaremba, Phys. Rev. A 68, 053611 (2003). 

[27] M. Modugno, C. Tozzo, and F. Dalfovo, Phys. Rev. A 70, 043625 (2004). 

[28] I. Danshita and S. Tsuchiya, Phys. Rev. A 75, 033612 (2007). 

[29] F. Dalfovo, S. Giorgini, L. P. Pitaevskii, and S. Stringari, Rev. Mod. Phys. 71, 463 (1999). 
[30] C. J. Pethic and H. Smith, Bose-Einstein Condensation in Dilute Gases (Cambridge University Press, 
Cambridge, 2002). 

[31] L. Pitaevskii and S. Stringari, Bose-Einstein Condensation (Oxford University Press, New York, 
2003). 

[32] A. Griffin, Phys. Rev. B 53, 9341 (1996). 

36 



[33] C. W. Gardiner and P. Zoller, Phys. Rev. A 58, 536 (1998). 

[34] R. Walser, J. Williams, J. Cooper, and M. Holland, Phys. Rev. A 59, 3878 (1999). 

[35] E. Zaremba, T. Nikuni, and A. Griffin, J. Low Temp. Phys. 116, 277 (1999). 

[36] H. T. C. Stoof, J. Low Temp. Phys. 114, 11 (1999). 

[37] H. T. C. Stoof, Phys. Rev. Lett. 78, 768 (1997). 

[38] N. P. Proukakis, K. Burnett, and H. T. C. Stoof, Phys. Rev. A 57, 1230 (1998). 

[39] H. Shi and A. Griffin, Phys. Rep. 304, 1 (1998). 

[40] L. P. Pitaevskii and S. Stringari, Phys. Lett. A 235, 398 (1997). 

[41] P. O. Fedichev, G. V. Shlyapnikov, and J. T. M. Walraven, Phys. Rev. Lett. 80, 2269 (1998). 
[42] S. Giorgini, Phys. Rev. A 57, 002949 (1998). 

[43] J. M. Cornwall, R. Jackiw, and E. Tomboulis, Phys. Rev. D 10, 2428 (1974). 

[44] J. S. Schwinger, J. Math. Phys. 2, 407 (1961). 

[45] L. V. Keldysh, Zh. Eksp. Teor. Fiz. 47, 1515 (1964). 

[46] P. Danielewicz, Ann. Phys. 152, 239 (1984). 

[47] J. Berges, e-print hep-ph/0409233 (2004). 

[48] T. Gasenzer, J. Berges, M. G. Schmidt, and M. Seco, Phys. Rev. A 72, 063604 (2005). 
[49] S. Konabe and T Nikuni, J. Phy. B: Atomic, Molecule and Optical Physics 39, S101 (2006). 
[50] K. Iigaya, S. Konabe, I. Danshita, and T. Nikuni, Phys. Rev. A 74, 053611 (2006). 
[51] S. Konabe and T. Nikuni, J. Low Temp. Phys. 148, 453 (2007). 

[52] H. Haug and A. -P. Jauho, Quantum Kinetics in Transport and Optics of Semiconductors (Springer, 
Berlin, 1996). 

[53] L. P. Kadanoff and P. C. Martin, Ann. Phys. 24, 419 (1963). 

[54] L. P. Kadanoff and G. Baym, Quantum Statistical Mechanics (Benjamin, N.Y., 1962). 

[55] J. E. Williams and A. Griffin, Phys. Rev. A 63, 023612 (2001). 

[56] J. E. Williams and A. Griffin, Phys. Rev. A 64, 013606 (2001). 

[57] R. A. Duine and H. T. C. Stoof, Phys. Rev. A 65, 013603 (2002). 

[58] M. Kramer, L. Pitaevskii, and S. Stringari, Phys. Rev. Lett. 88, 180404 (2002). 

[59] F. S. Cataliotti, S. Burger, C. Fort, P. Maddaloni, F. Minardi, A. Trombettoni, A. Smerzi, and M. 

Inguscio, Science 293, 843 (2001). 
[60] P. O. Fedichev and G. V. Shlyapnikov, Phys. Rev. A 58, 3146 (1998). 
[61] S. Giorgini, Phys. Rev. A 61, 063615 (2000). 

37 



[62] S. Tsuchiya and A. Griffin, Phys. Rev. A 72, 053621 (2005). 

[63] J. E. Williams, E. Zaremba, B. Jackson, T. Nikuni, and A. Griffin, Phys. Rev. Lett. 88, 070401 (2002). 

[64] S. Stringari, Phys. Rev. Lett. 77, 2360 (1996). 

[65] P. Navez and R. Graham, Phys. Rev. A 73, 043612 (2006). 



38 



